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Abstract 


Given  m  points  on  the  unit  sphere  in  n-space,  the  hyperplanes  tangent 
to  the  sphere  at  the  given  points  bound  a  convex  polytope  with  m  facets.  If 
the  points  are  chosen  independently  at  random  from  the  uniform  distribution 
on  the  sphere,  the  number  V  of  the  vertices  of  the  polytope  is  a  random 
variable.  We  obtain  an  integral  expression  for  EV  and  asymptotic  bounds 
of  the  form 


n  (n-6)/2,  , 

a  nv  "  (m-n) 


EV  S  Bnn(n-5)/2n, 
mn 
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I.  Introduction 


The  main  results  of  this  working  paper  are:  the  integral  expression 
(3.3)  for  the  expected  number  of  vertices  of  a  random  convex  polytope;  the 
upper  bound  (4.1)  and  consequent  asymptotic  upper  bound  (4.2)  on  the  expected 
number;  and  the  asymptotic  lower  bounds  (4.7)  and  (4.8). 

Briefly  and  informally,  the  bounds  are  these:  Let  , . . . be  inde¬ 
pendent  random  points  each  having  the  uniform  distribution  on  the  unit  sphere 
in  ,  and  let  V  be  the  number  of  vertices  of  the  (random)  polytope  whose 

facets  are  the  hyperplanes  tangent  to  the  sphere  at  . Then  there 

exist  constants  a  and  B  independent  of  m  and  n  such  that  asymptotically 

anntn-6>/2Cm-n)  s  EV^  s  Snn(n-5)/2m  . 

The  upper  bound  is  valid  for  any  n  >  2,  asymptotically  as  m  the  lower 
bound  is  valid  for  any  large  n,  asymptotically  as  in  +  ®. 

In  a  later  paper  we  will  present  the  results  of  numerical  approxima¬ 
tions  of  the  integral  (3.3)  for  moderate  values  of  m  and  n. 
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II.  Formulation 

T 

A.  Given  an  n-vector  a  and  a  real  number  b,  the  inequality  a_  x  <  b 

n  T  . 

defines  a  half -space  in  IRn  with  bounding  hyperplane  H  =  (x:  a/ x  =  b}.  The 

non-empty  intersection  of  a  finite  number,  say  k,  of  such  half-spaces 
defines  a  polytope  in  IRn  with  at  most  k  facets.  A  random  polytope  is  a 
polytope  generated  by  choosing  the  components  of  the  k  n-vectors,  a_p  • • • ,a^, 
and  the  k  real  numbers,  bp . . . ,b^,  from  some  probability  distribution  on  the 
real  line.  Obviously  the  structure  of  random  polytopes  will  depend  on  the 
choice  of  the  probability  distribution  from  which  the  data  are  drawn.  For 
examples  of  certain  distribution  schemes  the  reader  is  referred  to  [2,4]. 

In  this  paper  it  is  desired  to  consider  only  random  polytopes  in  Rn 
with  a  fixed  number,  m,  of  facets.  Hie  general  scheme  given  above  for  gener¬ 
ating  polytopes  will  lead  to  redundant  constraints  especially  in  the  case 
where  k  is  large  relative  to  n.  This  difficulty  can  be  avoided  if  bp...,b^ 
are  chosen  as  appropriate  functions  of  the  , . . . ,a^.  In  particular,  if  for 
j  =  l,...,m,  each  bj  satisfies 

(2.1)  b-  =  | a- |  =  (  X  a*)172  , 

J  J  i=i  J1 

then  the  random  polytope  generated  by  the  vectors  ap . . . ,  a^  is  circumscribed 
about  the  unit  sphere,  Sfi_p  and  has  exactly  m  facets.  Each  facet  lies  in 
one  of  the  hyperplanes  IT  =  (x:  a^x  =  |su  |)  and  is  tangent  to  at 

=  £u/  Iji..  | .  The  polytope  is  thus  completely  determined  by  the  unit  vectors 
£l , . . . ,  .  Consequently,  an  alternative  generating  scheme  for  this  type  of 

random  polytope  is  to  choose  the  vectors  from  some  distribution  on 

Sn-r 

Of  special  interest  here  will  be  the  case  where  are  independently 

and  uniformly  distributed  on  Sn_p  As  is  well-known  ([3])  this  choice  of  the 


II. 2 

tangent  points  is  equivalent  to  choosing  the  components  of  the  vectors  a^, 
j  =  l,...,m,  independently  from  the  standard  normal  distribution  on  the  real 
line  and  the  according  to  (2.1). 

It  should  be  pointed  out  that  requiring  the  polytopes  to  be  circumscribed 
on  the  unit  sphere  as  described  above  certainly  limits  their  possible  structure. 
Polytopes  of  certain  combinatorial  types  cannot  be  realized  with  all  of  their 
facets  tangent  to  a  sphere.  Therefore  it  may  be  argued  that  the  more  patholog¬ 
ical  polytopes  are  excluded  by  the  method  of  generation. 

B.  Let  Ei » • • • >Efn  be  independent  and  identically  distributed  points  on 

S  , .  As  discussed  above,  these  points  determine  a  unique  polytope  in  ]Rn 
n-1 

containing  Sn_i  and  having  exactly  m  facets.  Let  . . .  ,]^)  denote  the 

number  of  vertices  of  this  polytope.  The  function  Vmn  is  a  random  variable 
on  the  space  of  m-tuples  of  n-vectors,  with  values  in  the  nonnegative  integers. 
The  purpose  of  this  paper  is  to  investigate  the  expected  value  of  V^, 
specifically  when  the  have  the  uniform  distribution  on  S^. 


III.l 


III.  Integral  Expression  for  EVj^. 

A.  Fix  integers  m  and  n,  2  <  n  <  m,  and  let  £p...,j)m  be  independent 
random  points  chosen  from  some  distribution  with  density  g(£)  on  Sn_p 

Let  be  the  hyperplane  tangent  to  Sn-1  at  ap  i  =  l,...,n,  and  let 

P (£^ , . . . ,£^3  be  the  polytope  bounded  by  Hp . . . ,Hn«  Then  with  probability  1 

any  n  of  Hp.,.,1^  intersect  in  a  point  of  IRn  and  there  are  (™)  such  points 

of  intersection,  among  which  are  the  vertices  of  P .  ,2^) . 

Denote  by  A^  the  family  of  all  n-subsets  of  l,2,...^n.  For  any 

A  =  {i.,...,i  }  e  A  let  VA  be  the  event  that  the  point  of  intersection 
l  n  mn  A 

of  H.  is  a  vertex  of  P(£p . . . .j^) .  Then  V  is  the  sum  of  the  indica- 

1  n 

tor  functions  of  the  events  VA>  and  so 

Ev> '  aJ  Pr(VA’  • 

mn 

Because  the  £-  are  identically  distributed,  we  have  by  symmetry 

^nrn  '  ■ 

where  V  denotes  Vr.  . 

n  il,...,n> 

B.  Now  with  probability  1,  Ei>***»En  lie  on  a  unique  small  hypercircle 
on  Sn_1,  which  divides  Sn_1  into  two  caps.  Let  Cf^,...,^)  be  the  smaller 
of  these  two  caps.  Then  \’n  occurs  if  and  only  if  none  of  J^+p •  •  •  is  in 
C(Ei,..*,£n).  Given  Ei*,#*>En*  therefore,  the  conditional  probability  of 

Vn  is 

PrtVn|£l . £„)  -  (1  -  /  g(E)d£)m'n  . 

C(Ep  •  •  •  *En^ 

So  EV^  is  (™)  times  the  expected  value  of  this  function  of  Ei» •••»£«*  t*iat  *s» 

^mn  =  O^^nlEl . =  _  /  **•  /  ^  ‘  f  /  steW"1'" 

Sn-1  Sn-1  CvEl»*”»2n) 

•  d^,...,^  .  (3.1) 
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C.  Now  we  assume  that  the  common  distribution  of  i-s  the 

uniform  distribution  on  S„  In  this  case 

n-i 

area  of  C^,...,]^)  In_2M 

g(£.)d£  area  of1  S  I  21  -,(tt/2)  * 

n-1  n-2 

where  r  =  rC^,...,^)  is  the  angular  radius  of  (0  <  r  <  m/2) 

and  where 


I.(r)  =  /  sin*Sc  dx  (0  <  r  <  tt/2,  k  -  0,1,2,...) 

0 

is  the  area  of  a  cap  of  angular  radius  r  on  S^.  Thus  (3.1)  can  be 
rewritten 


EV  =  (®)E 
mn  n 


1  - 


21 


In-2Cr)  fn 

W7)\ 


n-2 


where  the  random  variable  r  is  the  radius  of  C(£^, . . . . 

Now  it  follows  from  results  of  R.A.  Miles  ([5],  theorem  4,  p.  368) 

2  2 

that  t  =  sin  r  has  the  beta  ((n-1)  /2,  1/2)  distribution,  whose  density  is 

2 

proportional  to  t^n-1-*  ^2^’1(l-t)"1//2  on  (0,1).  From  this  it  follows  that 
the  density  of  r  is  fn(n_2)(r)»  w^ere 

fk(r)  =  C^sin^r  ,  0  <  r  <  tt/2 


and 


°k 


2r(^) 

r(y)rc^j4 


k  =  0,1,2,...  . 


(3.2) 


Moreover ,  the  distribution  function  of  this  density  is 

Ik(r) 

Fk(r)  ■  -  0  •= r < ,/2  • 


Using  the  above  notation  we  obtain  the  expression 


III. 3 


_  tt/2  . 

^*■0/  a-K-2M>m"\c„-2)(r)dr- 


(3.3) 


Although  we  will  not  use  it,  the  following  alternative  expression  may  be  of 
interest.  If  and  denote  the  distribution  function  and  density  of  the 
beta(k/2,  1/2)  distribution,  then  the  substitution  t  =  sin  r  in  (3.3)  gives 

“■->  '  0  p  -  7  Vi<t»"'"Vi)2Mdt  • 


IV.  Upper  and  Lower  Bounds  jy  ^ 

A.  The  evaluation  of  (3.3)  is  not  difficult  in  case  n  is  2  or  3: 


EV 


^  V2„  r.m-2  2 


m2 


=  (?)  /'  (1  -  ^)m"^  —  dr  =  m(l  -  (£)m_1)  , 

=  (?)  /  (1  -  ^-(1  -  cos  r))m"  •  j  sin3r  dr 
=  2m- 4  -  (j)m-1(m+l) (m-2)  . 


and 


EV 


m3 


(We  note  in  passing  that  m  and  m-1  are  the  only  possible  values  of  Vm2, 
the  latter  occurring  only  when  » •  •  •  are  all  in  one  semicircle.  The 
value  of  Em2V  given  above  gives  the  well-known  value  m(-^)m  1  for  the 
probability  of  this  event.) 

The  above  calculations  suggest  that  for  fixed  n,  EVmn  is  asymptotically 
of  the  form  Km.  This  section  is  devoted  to  finding  upper  and  lower  bounds 
on  EV^  which  establish  this  asymptotic  growth  rate  and  bound  the  growth 
of  for  large  n.  (Unfortunately,  the  tempting  conjecture  that  =  n-1 
is  false:  is  approximately  of  order  (Bn)^2  as  we  have  noted.) 


where 


B.  First  we  will  establish  an  upper  bound  by  showing  that  for  2  <_  n  <  m, 

“mn  s  -  V  W-D 

>n-3 


a  =  c  r-L-^lndr 

Ln(n-2)  LCn_2J 


mn 


Before  proceeding,  we  remark  that  B  is  negligible: 


Bnm  <  • 


-£f  <  T  when  n  2  3>  %  s  C^)1"'1. 


I 


and  since 
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We  observe  also  that  Wallis's  formula  ([1],  p.  258),  viz. 

r(a+  h 


implies 


and 


?  1„  1  .  1^ 
frS7TT~~(1  “  H  +  0(t3)  as  a 


& 


Cn(n-2)  ~  ^  n 


a 


~  /2tt/  [n"2)  (1  ‘  Tfri-T)'-*  * 


Thus 


.  ..  (2n)n/2  (n-5)/2 

\  ;;  T/T-  n 

Tie 


and  so  if  m  £  n,  then  asymptotically  as  n  +  « 


The  proof  of  (4.1)  begins  with  (3.3),  which  says  that 


^mn  0("n(n-2)^n-2,m-n,n(n-2) 


where 


tt/2 


sj u  =  [  a-7Fj(r))ksinJlrdr- 

If  l  >  j+1,  an  integration  by  parts  using 

u  =  _  jr-  sin^  -*r  ,  dv  =  (1  -  ^Fjfr))^-  -j  CjSin^r)dr 

*/2 


yields 


(4.2) 


(4.3) 


(4.4) 


*,*  *  JO  -  cos  r  dr 


and  hence 
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2  1  l.k+1  2  &-j  n 

C7FTC7J  C7  Tc+T  Sj  ,k+l,£-j-l 


for  £>j+l  . 


(4.5) 


Iterating  (4.5)  gives,  for  Z  >  v(j+l). 


Sjk£ 


-l)j-(i-2))  (-l-jk+i 


r  2/(Jl-j)(^2j-l) 

T7J  — Tk'-fTTOc^' 


(il-vj-(v-l))  s 

Tk+vT  j  ,k+v,J.-v(j+l) 

(4.6) 


In  particular  this  is  valid  for  j  =  n-2,  k  =  m-n,  Z  =  n(n-2),  and  v  =  n-2. 
In  this  case  for  i  =  2, 3,..., n-2  we  have 


(£-j)  (8.-2j-l)  . . .  (£-(i-l) j-(i-2)) 


=  [  (n-1)  (n-2)  ]  [  (n-1)  (n-3)  ] . . .  [  (n-1)  (n-i)  ]  =  (n-l)1"1 
and  also  by  a  simple  integration 


S, 


=  S 


2  1 


’j  ,k+v,£-v(j+l)  n-2,m-2,n-2  Cn_2 


Putting  the  two  expressions  above  into  (4.6)  yields 

S  <  f_  2  I11"1  (n-2) !  (m-n) !  f  ,  ->n-2 

n- 2, m-n, n (n-2)  ^n_2J  (m-1) !  ^ 


lsm-1^ 


rr(l  -  (4T  A) 


ny1r  2  -i  (n-2)! (m-n)! 
i=ln-2J  (n-i-1) ! (m-n+i) ! 


(n-l)1"1 


,l^m-n+i 

¥ 


Ntiltiplying  by  Cn^n_2^  (Jj)  gives  an  upper  bound  on  EV^;  simplifying  and 
letting  the  sum  run  over  j  =  n-l-i  then  gives  (4.1). 


C.  Next  we  show  that  for  any  n  >  2  and  any  s  e  (0,1), 

EV  >.  -fo  1--'  - — (n-l)n  ^(m-n)  ,  asymptotically  as  m  +  ®  . 

m  dn",  nle11"1 
n-l 


(4.7) 
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where 


n  /'i  + 


d.  = 


3  r(£>r(^) 

Before  proceeding  we  notice  that  Wallis's  formula  implies 


dj  ~  /2/nj  as  j  -*  <*>  . 


This  and  Stirling's  formula  provide  that 

n-l  _  c 
n-o 

2 

e  n  as  n  +  “  . 


d(n-l)2  2n~1en  ^n_1-jn-l  _  (2tQ  pn  n  2 


,n-l  ,  n-l 

dn-l  n!e 


ite 


i rr 


Therefore,  for  fixed  large  n,  (4.7)  implies  that  for  any  e  e  (0,1),  we 


have  asymptotically  as  m  +  « 

n-l 


n-6 


EmnV  -  ~2~T7?  ^  n  2  *  (4-8) 

Tre 

To  prove  (4.7)  we  begin  again  with  (3.3).  For  any  a  e  (0,y) 

f<1-TFn-2t“»'"'nf„(n-2)Wdr 
■  0<1-7F„-2M)""Fn(,-2)W  •  W-9) 

It  can  be  checked  that  for  0  <  a  <  j 

1  lc*f* 

Fj, (a)  =  cos  a(dk+1sin  a  +  dk+;jsin  a  +...) 

by  differentiating  and  using  elementary  properties  of  the  gamma  function. 
Moreover,  dQ  >  dk  >  d2>...  .  Therefore 


(dn-l(sinn"la  +  sinn+la  +•••)) 


1 


cos  a 


and 


IV,  5 

d  ..sin11 

1  -  -2zk - 

i.  cos  a 


Fn(n-2)Ca)  "  d(n_1)2  cos  a  sin(n_1)  a  . 


So,  writing  t  for  sin  a,  we  have 


4  0d,  .,2  ^  ** 

(n-1) 


-l  - 

n-1  t 


n-l'im-n 
==  for  any  t  e  (0,1). 

-T2' 


/1-r 


Let  e  be  an  arbitrary  number  in  (0,1);  then  we  have 

2  0d  n2  ^tn'l)^l  -  %  A1/-  1£  o  <  t  <  j 

(n-1) 

and  0  <  x  <  /1-c2 

d  i  i 

If  we  now  write  a  for  —r~—  t  ,  then  this  inequality  becomes 


V  2  0dr  112<r2->n'1  if  0  <  ,  <  1 

(n-1)  n-1 


n-1 


and  0  <  a  <  -^lifl-e2)  4  , 


Notice  that  the  value  of  a  maximizing  the  right  side  of  the  above  inequality 


•  _  n-1 
m-1 


is  ~  .  So  if  24.  <  4£kl-eVn"1)/2 


,  then 

c  \r  a(n-l)2  0n-l_n,n-Un-l,.  n-Km-n 

bmnV  -  V  fa - 2  e  fa)  (1  '  fa 

n-1 

Since  (JJ)  >  , 


m-l  lz 

i< 

dr 


Em,V  2  (m‘n)^ri-  ^  Cl  -  , 

n-1 


and  (4.7)  follows  from  this. 
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points  are  chosen  independently  at  random  from  the  uniform  distribution  on  the 
sphere,  the  number  V__  of  the  vertices  of  the  polytope  is  a  random  variable. 
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We  obtain  an  integral  express ion vf or  EV_>  and  asymptotic  bounds  of  the  form 
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